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ABSTRACT 


This study demonstrates the significance of water vapor's 
influence on sSatellite-retrieved aerosol characteristics using 
NOAA's AVHRR. An improvement to optical depth and Aerosol 
Particle Size Index (S12) estimations derived from channels 1 
(0.634m) and 2 (0.86um), is made through knowledge of column 
water vapor derived from channels 4 (10.8um) and 5 (12.0um). 
A 2.0 gm cm-2 column water vapor produces a 5% increase in 
S12. This results in a 15% error in the variable scattering phase 
function P(@) and retrieved aerosol optical depth (6g). The error 
introduced by water vapor is quantified through use of the 
LOWTRAN7 atmospheric propagation model to be applied as a 
reformulated parameterization of (P(®)). 


il 


I. INTRODUCTION. ........0:...ccccccnncccssscecsosecinmaceessccoccssesmnctaam l 
A_ PURPOSE ooic.....ccccccceoeoocescotteceiiiaiaincsssocscasenencresesseseceseasec are slacmiuaa l 
B. MOTIVATION. ....cccecccccccosccccnscascccescncssescsscnttencnsccsetensracrsesacs simi 2, 
C. THE MODEL. .......0.:000ccescccserentcneesesscnaeerecencsccnssovecnvsocosteestaeenn iia 3 
TT. THEORY cccccoccecccttiies sss cccccsineeeseeslersmeciiesiss™= re ensere> ssn sl ams S 
A. PARTICLE SIZE INDEX (S12)........:sccseccneeceeeereeeeeseneceeccnsersscneeeees 5 
B. WATER VAPOR ABSORPTION IN THE VISIBLE / NEAR 
INFR A-RED......0:cccecccceoosessatee «stein seccueees-auecuumnmemenegar++-- 00+" ¥camau 10 
C. BRIGHTNESS TEMPERATURE DIFFERENCE (14-15).........eem 14 
LIL. PROCEDURE. z.......ccccccccccccccvccssccccccccccccssssccccsssscssscccseesess 20 
A. INTRODUCTION .......00-.00sccrecccssssscseeccaseeeessssceaemanernyesccss= eras deaaam 20 
B. OVERVIEW... ccccccsescesssccecceteetcesssccsscncscseeecensssnenminccronssss eases acugiam 20 
C. INPUT CHARACTERISTIGCS...........:::cecccecsseseeeceeneentessesseenaneeser sees 24 
1. Optical Depth Profiles...............sssscccssercerssseseesenneneneerersecesnens 2 
2. Water Vapor .........0:--cssesscscecnrececocescesvenerssenecennetacnsencrevic ssa 26 
3, Surface Temperature...........cccceceeeessssessrtenneeesceseeseesersesesses eens 28 
D. PROCEDURAL SUMMARY ........::csescccseecneseenseensecueceuneeeanenaeeren ees 28 
TV, RESULTS cccccsccecccctecscseusdeesestess+sss-ccrsse+s => sss ces smneenocces maa 30 
A. THE EFFECT OF AEROSOL OPTICAL DEPTH VARIATIONS 
ON S12. oceccecocecdecsonecscthesteneossne sees tsumsanteessetciiaesre-ccmeennnocsr eam 30 
1. Multiple 5g(A) Scheme ............ssscseserseeseecesseesesrarentneestnesssrens 30 
2. Parallel Sa(A) SCHEME ............cscccceceescsesssessseeeseeeteeeeecsetsntee aeons 30 
B. THE EFFECT OF SURFACE TEMPERATURE VARIATIONS ON 
THE T4-Ts5 REPRESENTATION OF W........:seseseseeseessrseeeeecetesress 30 
C. THE W CORRECTION TO S812......::::: eeeeeseeseseesereseseneceeneeeeesees ones 34 
1. Graphical Application ............::ccssscsessseeectessseeensneesttneres testes ees 34 
2. The Reformulated Parameterization. ............ssseeerereseseeseeerseseess 36 


TABLE OF CONTENTS 


De -SINGCE VERSUS MULTIPLE SCAT TERING 22.205. ch eae ooo 


V. CONCLUSIONS 


@ee@eeeeecoeeegc ee eee ee eee eeeeeeze eo eeeeeeeeseeeoeeeeeeeeeoseenceveeneueeeeeee eee 


LIST OF TABLES 


Table 1 NOAA-7 AVHRR Channel Bandwidths................. cc cccccccceescceceeceeeeeess 5 


Table 2 Model Output Coefficients Applied To Dalu's Theory 
(CW SACTA-15)).....cs.000 025 e-aagenwcteaeec tenets ssiesan ei eee ee nee ocr 34 


vl 


LIST OF FIGURES 


| Extinction formulation as a function of particle radius for channels | 


(solid) and 2 (dash). Modified from Durkee (1984). 000.0... cecececeeeees 9 
.2 NOAA-7 AVHRR channel | and 2 response functions and total 

‘SGqOVSSS 0) 1 ley g (ei #: (0.92 (el oc ns on a oe ene ene ne Ree eee 1] 
. 3 NOAA-7 AVHRR channel | and 2 response functions and atmospheric 

MAliGiMitancerGUc 1O-Waleh Vv ADOF.«.es.cocs.420. cine. 2 utaasscne oes se ecneeeae, vacen ae 12 
.4 NOAA-7 AVHRR channel | and 2 response functions and the 

influence of water vapor on atmospheric radiance. ..............cceeeeeeeeeeeeeees i 
.5 NOAA-7 AVHRR channel 4 and 5 response functions and atmospheric 

Uhalistmiliteance<cdue tO Waler VANOL..:. 22.25. -2acseasee- ncn .-seees coese ae ee rcien ea aee 16 
.6 Total water vapor content (kg/m2) as a function of brightness 

remperature diterence trom Daly (1986)... jocic-ecsincisstausteseseceese cose ss cau 18 
. 7 Comparison of water vapor content from ships with radiative transfer 

model runs using ship radiosonde data. From Dalu (1986). ................... ie, 
MOM ce Mat Ol AMOMe MMM ac aaiancicdeeaie ieee eens es oak eee ee ect 21 
.9 NOAA-7 AVHRR filter function reSPONSE CUIVES............ ce eeeeceeeeee eee e es ZS 
Oma Uitipie Optical dept protiles (VE) ic .cr..eccc-ccsnen swoon sousecesincssoadseoasens 26 
meme iearallelkoptical-depth- proiiies (DP) 20.5 ner synscssssccceeceeue see owen settee se asee 27 
. 12 S12 versus column water vapor (w) for the multiple optical depth 

OMENS ACNAN) rrctesetea cares nie seine cece neeiniie waite tase stta see natee ecee Se Oana oe Cums « caleae reanmie ees 31 
. 13 S12 versus column water vapor (w) for the parallel optical depth 

SENET 0 Ce ene a ee nn en er Seen emer corre etree Tr er nee 32 
. 14 T4-T5 versus column water vapor (w) for a variety of surface 

BOM ORACLE Socata ste he eee eee cue one ese een aceon nn cne ethene lett macehee tree 33 
mrs 51) COMECHONMOL COMIMM Wale VAPOR... 612s... .. cosy vereeee essence sere cee ees 35 

om onPanOrmializea DY. Ol Ky) eae creas sas nnsnecnnegscsessc sae-e ss ecngeesneuee areas 37 


ACKNOWLEDGMENTS 


The author wishes to recognize the key players, without whom, this endeavor 
might have become a struggle. The top of this list has to be my wife, Mary, who 
has dealt with many of the worlds problems, so that I didn't have to. I am 
grateful to: Professor Philip A. Durkee, for his consistent demeanor and 
patience despite my continual knocking at his door; Professor Kenneth L. 
Davidson, for particularly quick turnarounds as a reader; and Professor Curtis 
A. Collins of the Oceanography Department, for allowing me to use his 
computer. I want to be sure to thank Mr. Craig Motell, for always being 


available for programming support. 


Vill 


I. INTRODUCTION 


A. PURPOSE 

Direct measurements of aerosol characteristics can be provided by ship or 
aircraft. Although they provide valuable information, they are very expensive, 
require detailed coordination, and only small sections of the world's oceans can 
be monitored at any time. Satellites can provide near continuous spatial and 
temporal data desirable for evaluating large scale global effects as well as timely 
local measurements suitable for military applications. Meteorological satellites, 
however, have been designed to measure cloud albedos while much lower 
radiances are associated with aerosol optical depth variations (Durkee 1984). 
Therefore, careful processing of satellite data is necessary to determine aerosol 
characteristics. 

Frost (1988) has developed an algorithm to derive aerosol particle 
characteristics from the Advanced Very High Resolution Radiometer (AVHRR) 
on the NOAA polar orbiting satellites. By taking advantage of the varying 
Sensitivities of aerosol scattering in two spectral windows, a multi-channel 
comparison provides information about the aerosol particle size distmbution. By 
applying this information to a variable scattering phase function, improved 
estimations of optical depth are possible. 

This study will provide an improvement to Frost's technique. Frost requires 
analyzed pixels to be free of clouds. However, these clear pixels still contain 
varying amounts of water vapor. Frost assumes that varying amounts of water 


vapor do not appreciably effect retrieved radiances. This study will use the U. S. 


Air Force's LOWTRAN7 atmospheric propagation model to demonstrate the 
degree that water vapor influences present aerosol estimations. 


This thesis will: 
1. Show the effect of varying water vapor amounts on satellite- 
retrieved aerosol characteristics, specifically the particle size index 
parameter or S12 as defined by Frost (1988). 


2. Provide a parameterization of S42 corrected for column 
water vapor (WwW), using a multichannel analysis of AVHRR data. 


B. MOTIVATION 

Radiation in the atmosphere is attenuated by many processes. It is absorbed 
and scattered by molecules, aerosols, and meteorological features such as clouds, 
fog, and precipitation. Some of these influences may be of short duration and of 
local scale such as a cumulus cloud passing overhead and casting a shadow. 
Others, such as stratospheric dust particles are semipermanent and of global 
proportions. 

Aerosols are solid or liquid particles dispersed in the atmosphere. Dust, 
smoke, and even fog or rain are examples of aerosols. Aerosol particles are 
present in the atmosphere in various forms and concentrations. In this thesis, the 
change in particle size due to coalescence of water vapor will not be considered. 

With the ever increasing use of electro-optical devices, military applications 
benefit significantly from improved aerosol measurements. The optimum 
waveband for a laser is from 0.4 to 1.0 um. This waveband avoids ozone 
absorption at shorter wavelengths, as well as the water vapor absorption that 
dominates longer wavelengths (Bloembergen, et al 1987). This waveband, 


however, is significantly influenced by aerosol scattering. 


Another important hypothesis that demonstrates the need for better aerosol 
estimations and modeling on a much larger scale is that of Charlson, et al. 
(1987). Charlson proposes that the major source of cloud condensation nuclei 
(CCN) over remote unpolluted oceans appears to be aqueous dimethylsulphide 
gas (DMS) produced by phytoplankton as a waste product. DMS is emitted 
through the sea surface through a chain of complex biological and chemical 
processes where it is oxidized in the atmosphere to form sulfate aerosol particles 
which act as CCN. Increased phytoplankton populations and the resulting 
increase in aerosol particles effectively cool the atmosphere by increasing 
backscatter to space. A 30% increase in aerosols over the ocean surfaces 
worldwide can conservatively decrease the average global temperature by 1.3 K 
(Charlson 1987). 

For many applications, aerosols result in noise interfering with 
measurements of other processes. These applications benefit from improved 


aerosol estimations that eliminate aerosol contamination of the data. 


C. THE MODEL 

LOWTRAN7 (Low Resolution Transmittance) is a model that calculates 
atmospheric transmittance, scattering and absorption. It incorporates separate 
molecular profiles for all major as well as 13 minor and trace gases, all as a 
function of altitude (0 to 100km). Representative atmospheric aerosol, cloud, 
and rain models are provided with the ability to replace any or all values with 
theoretical or measured data at the discretion of the user. The model also 
provides six reference atmospheres with options to redefine or retain profiles. 

LOWTRAN7 is a band model. Average transmittance is calculated over 20 


cm-! bands. The band parameters have been determined semi-empirically 


determined semi-empirically through the averaging of line parameters. In 
principle, the line by line method is available, but the calculations are too 
elaborate and time consuming for most users. Some accuracy is gained by a line 
method, but it is not significant enough to warrant the loss of efficiency. 
Individual lines within a band may vary significantly, yet it is fair to characterize 
the effective band parameters through the averaging of line parameters so long 
as a specific line result is not desired. 

For use in multiwavelength applications (ie. satellites) the appropriate 
spectral response functions for the sensor can be applied to the output data 
through LOWFIL7 (LOWTRAN?7's filter program). The data set derived 
through the running of the main program is filtered through the sensor's 
response function and produces a normalized radiance estimate of the sensor 


measurement. 


Il. THEORY 


A. PARTICLE SIZE INDEX (S12) 
The Advanced Very High Resolution Radiometer (AVHRR) has been a part 
of all the NOAA polar orbiting satellites. It has five channels: 


TABLE 1 

NOAA-7 AVHRR CHANNEL BANDWIDTHS 
Channel # Bandwidth (im) 

I 0.58 - 0.68 

2 0.725 - 1.10 

3 S55 9.05 

4 10.30 - 11.30 

5 11.50 - 12.50 


AVHRR satellite detected radiances in the red-visible (chan. 1) and the near 
infra-red (chan. 2) wavelength ranges have been found to be positively 
correlated with aerosol optical depth (6) (Durkee 1984; Griggs 1979). 

Using a single scattering approximation and knowing that optical depths over 
the ocean are small (Pfeil 1986), Frost (1988) used a simplified form of the 
Radiative Transfer Equation (RTE): 





L(3,n,9) ~ Po? p(@)3 (2.1) 


where, L = diffuse radiance, 


Wo = single scatter albedo, 


Fo = solar radiative flux, 
ut = cos@ (6 = satellite zenith angle), 
p = Scattering phase function, 
© = single scattering angle, 
5 = optical depth. 
In an optically thin (small 5) atmosphere, multiple scattering effects are very 
small, making single scattering a reasonable assumption (Frost 1988). 
There are three primary sources of upwelled radiance (L) that contribute to 


satellite detected radiance in the red-visible and near-infrared wavelengths. The 


contributions can be expressed as: 


L=La+Lr+Ls V5) 
where La is aerosol scatter, Lr is molecular or Rayleigh scatter, the wavelength 
of radiation is much greater than the size of the scattering material, and Lg is sea 
surface reflectance. Except for sunglint, Ls is very small for red wavelengths 
(albedo = 0.5%) and zero for wavelengths greater than 0.7m (Ramsey 1968). 
Therefore, as long as sunglint geometry is avoided, total radiance is due to the 
sum of aerosol and molecular scattering effects (L = Lg + Lr). 

Lr is strongly wavelength (A) dependant, but does not vary spatially. 
Therefore Lr forms a baseline L(A) and aerosol particles are the dominant 
source of upwelled radiance variations (AL ~ ALa). From Eqn. 2.1, changes 
in satellite measured radiance are directly determined from changes in aerosol 
characteristics like Wo, p(©), and aerosol optical depth (6a). Therefore, to 
retrieve a more accurate 5 from Lg measurements, knowledge of p(®) 


variations 1s required. 


Durkee (1984) showed that retrieved radiances at red wavelengths (chan. 1) 
are more sensitive to varying densities of small aerosol particles than that of 
near-infrared (nir) wavelengths (chan. 2). The most significant parameter 
contributing to retrieved radiances at AVHRR channel | and 2 wavelengths is 
the extinction coefficient (Gex¢) and specifically the vertical integral of Oext or 


the optical depth (6). 


H 
h= Oext dz (2.3) 


Oext 1S the sum of extinction due to absorption (Oabs) and due to scattering 
(Oscat). For marine aerosols, absorption is negligible, therefore Gext = Oscat 


and 


H 
o= Oscat dz. (2.4) 
Three terms contribute to Oscat: 
Oscat = 4 rr? Qscat(m.d.r)m() dr vs) 


where, r= particle cross sectional area, 
Qscat(m,A,r) = scattering efficiency, 
m = complex index of refraction, 


dX. = wavelength, 


r = particle radius, 


n(r) = number distnibution of particles. 


Figure | shows the influence of the terms of Eqn. 2.5 as a function of r and 
the resulting number-weighted extinction (mr2-Qscat-n(r)) for channel 1 and 
channel 2. Scattering cross section (mr2-Qscat) and the particle size distribution 
(n(r)) for a modeled distribution of marine particles at 80% RH are plotted. As 
wavelength increases, so does the r of most efficient scattering (Qscat). The 
result is a shift in the scattering cross section and a similar shift in the r that most 
significantly contributes to Oscat, and therefore L at a given A. Therefore, the 
slope of L(A) closely resembles the slope of 5(A), and is similar to the slope of 
the aerosol particle size distribution (n(r)). Since aerosol particles most 
effectively scatter radiation at a wavelength similar to the radius of the particle, 
an increase in the number of small aerosol particles will increase the ratio of 


Lehi to Lch2, and the spectral variation can be quantified: 





Wo: F 
S—*.9(®)-5(ch1) 


L(red) _ Lychl) _ _4 _ (red) (2.6) 
Linir) Lich2) ®o:Fo O(nir) 


~ gy P(O)-0(ch2) 


This ratio is called the Particle Size Index (S12) (Frost 1988):! 





Ccatinil 


_ (D)red (O)red 
512= Dir ~ Ohnir oa 


\While Frost removed Lr prior to taking the ratio, it will not be removed in 
this study. It shall be accounted for in the parameterization of p(®). 
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Fig. 1 





0.001 


0.1 10 
Rodius (um) 


Extinction formulation as a function of particle radius for 


channels 1 (solid) and 2 (dash). Modified from Durkee (1984). 


Satellite measured radiances therefore cannot provide the actual particle size 
distribution, but the S42 gives the slope of the size distribution curve. In 
general, the greater the contribution of smaller particles, the larger the ratio. 

Figure 2 shows a plot of NOAA-7 AVHRR channel | and channel 2 
response function curves superimposed on LOWTRAN7's Tropical model 
atmosphere output radiances (L(A)). Although as this paper will show, water 
vapor absorption does effect S12, the dominant variable influencing L(A) in 
these wavelength bands is aerosol scattering. For the red and near-infrared 
region of the spectrum, variances in L(A) are proportional to changes in da(A) 
and therefore aerosol scattering. Note that L(A) in the channel 1 window (short 
wavelengths, small particle scattering) is greater than that of channel 2. This is 
primarily due to the slope of Ly(A) and leads to an Sj? ratio that is generally 


greater than 1. 


B. WATER VAPOR ABSORPTION IN THE VISIBLE / NEAR 
INFRA-RED 


Although aerosol scattering is the dominant influence on satellite retrieved 
radiances at visible and near-infrared wavelengths, the role of atmospheric water 
vapor is not insignificant. The magnitudes of the overall L in these wavelengths 
are not appreciably effected, but often the goal of using satellites to estimate 
aerosol particle characteristics is to look for variations, which are an order of 
magnitude smaller than the total L. 

In order to observe sometimes subtle, yet often significant trends such as 
global climate variations, an improvement in the accuracy of aerosol estimations 


is necessary. 
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Fig. 2 NOAA-7 AVHRR channel 1 and 2 response functions and total 
atmospheric radiance. 


Figure 3 shows the effect of water vapor alone on t for AVHRR channels 1 
and 2. Notice that for channel 1, t is greater than channel 2, meaning that less 
absorption is taking place. This would tend to decrease L2 more than Lj and 
thus lead to an over-estimate of the S32. Therefore a change in water vapor in a 


region may be misconstrued as a change in aerosol characteristics or disguise a 


significant change. This implies that we might expect to see a minimum S}12(w) 


at W = zero, and an increase in S12(w) with w. 
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Fig. 3 NOAA-7 AVHRR channel 1 and 2 response functions and 
atmospheric transmittance due to water vapor. 


In Figure 4, the influence of water vapor is again demonstrated when water 
vapor is added to LOWTRAN7's Tropical atmosphere model. Notice that 
while radiances decrease in both channel wavebands, channel 2 is more 


significantly effected. This again will cause S12 to increase with w even though 
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Ly and L2 decrease. There does appear to be some increase in L below 0.6 Lm. 
This is most likely due to an effective shift in particle size due to coalescence of 
water vapor. This study intends to look at the influence of water vapor alone on 


retrieved radiances, and aerosol growth will be inhibited. 
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Fig. 4 NOAA-7 AVHRR channel 1 and 2 response functions and the 
influence of water vapor on atmospheric radiance. 
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C. BRIGHTNESS TEMPERATURE DIFFERENCE (Tq4-Ts) 

In order to take advantage of the water vapor parameterization of the 
scattering phase function proposed in this study, some method of measuring the 
water vapor present would be required. 

A significant benefit of using the AVHRR is that while channels | and 2 
provide information on optical depth variations, channels 4 and 5 in the 10 to 13 
l1m window region have demonstrated usefulness in atmospheric water vapor 
retrieval. 

Unlike the aerosol backscattering of channels 1 and 2, the brightness 
temperature is the result of single direction radiation transfer, from the sea 
surface to the sensor. Therefore, Tt is directly proportional to retrieved T, with 
absorption reducing the radiance reaching the sensor. 


From Dalu, et al. (1981): 


Ts = Ta + g(w)(Ta-TB) (2.8) 


Ts = Sea Surface Temperature (SST) 
Tq = Brightness Temp. of most transparent channel 
Tp = Brightness Temp. of second channel 
g(w) = function of water vapor content 
1 - Te 
& (ea ai) 2 a 
where Tq and tf are total atmospheric transmittances for channels o and B and 


C=(Ts- T B)(Ts- Tq). Solving for (To-Tp): 


=| 


(Ts-Tq) 
(Ta-Tp) = an (2.9) 
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By using an approximation proposed by Prabhakara, et al.(1974): 
(Ts-Tq) = (Ts- T g)Kw secé (2.10) 


where @ is the scanning angle and K is an absorption coefficient. Column water 


vapor (w = Jpwdz) (g m-2) can be related to the brightness temperature 
difference: 
(Ts- Ta) 


(Tq-Tp) = K a(w) w sec@ (2.11) 


Although both g(w) and (Ts- a) increase with increased water vapor, their 


ratio remains nearly constant (Dalu 1986). This leaves us with: 
w = A(Tq-TB)cos® (2.12) 


where A = const. = (g(w)/(K(Ts- T q)). 

To maximize the effectiveness of this scheme, Tq should be minimally 
influenced by water vapor in the atmosphere to allow the retrieval of water 
vapor temperatures near the surface. TI, on the other hand, should be more 
sensitive to the atmospheric water vapor distribution. Under such conditions, 
(Tq-TB) would indicate the influence of water vapor, which would be 
proportional to the amount (w). 

Channel 4 (~11um) and channel 5 (~12~m) windows of the AVHRR were 
designed for sensing of the sea surface temperature (SST) corrected for 
atmospheric absorption (Dalu 1986). Despite being optimized for SST retrieval, 
they describe a linear function from which the atmospheric water vapor 


correction can be derived (McMillin and Crosby 1984). Figure 5 shows 


NOAA-7 AVHRR channels 4 and 5 response function curves supemmposed on 
three model atmosphere t of varying water vapor amounts. Notice that Tq is 


generally greater than 15, leading to a positive T4-T5, and that higher T4-T5 


correlates directly with higher water vapor densities (Tropical > US Standard > 


Midlatitude Winter). 
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Fig. 5 NOAA-7 AVHRR channel 4 and 5 response functions and 
atmospheric transmittance due to water vapor. 
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Although not the ideal sensor windows for water vapor retrieval (the 


channels are too similar), the use of T4-T5 is justified to within the accuracy of 
the approximations made to achieve Eqn. 2.12. The most beneficial aspect of 
T'4-Ts5 is the co-location and temporal linkage with channels | and 2, providing a 
water vapor snapshot simultaneously and using the same path as the aerosol 
optical depth (6) measurements. 

Therefore, an estimate of vertically integrated water vapor density (w) is 


provide by: 
w = A*(T4-T5)cos@ (2.13) 


Using A* = 19600 g /(°K m2) and a scanning angle 0 = 0, Dalu (1986) applied 
this relationship to a radiative transfer model using a variety of temperature and 
humidity profiles. Figure 6 shows the correlation coefficient for the resulting 
line is R = 0.99 with the error given as + 1.5 kg m-2. Comparing results 
calculated using radiosonde data with shipboard measured water vapor amounts, 
Figure 7 shows the correlation coefficient dropped to R = 0.78 with an error of 
+ 40kg m2. The decrease in accuracy most likely resulting from calibration 
errors as well as inexact temporal and spatial correlation of ship and satellite 
data. 

At channel 4 and 5 wavelengths, scattering is a relatively insignificant 
influence determining brightness temperature (T). The dominant factor is water 


vapor absorption of surface temperatures. 
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Fig. 6 Total water vapor content (kg/m2) as a function of brightness 
temperature difference from Dalu (1986). 


Some of the long wave radiation emitted by the sea surface is absorbed by 
water vapor. The vapor then re-emits radiation based on temperature in 
accordance with the Planck function. If the atmosphere, and the water vapor 
included, were not generally cooler than the sea surface, the absorbed radiation 
would be re-emitted at the same rate, not having any effect on radiance or 


brightness temperature. 
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Fig. 7 Comparison of water vapor content from ships with radiative 
transfer model runs using ship radiosonde data. From Dalu (1986). 
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IIT. PROCEDURE 


A. INTRODUCTION 

This study uses the LOWTRAN7 model to improve algorithms used in 
satellite retrieved optical depth measurements, specifically a correction due to 
varying amounts of water vapor present. The benefit in using a theoretical 
model is that the parameters involved can be specifically quantified and 
controlled. Therefore, a precise relationship between satellite measured 
radiances (L, w/cm2-steradian-cm~!) and varied parameters such as optical depth 
(d), atmospheric temperature (T, °K), column water vapor (w, gm m-2) can be 
determined. These relationships will be consistent with that of the real 


atmosphere to the degree that the model represents the real atmosphere. 


B. OVERVIEW 

Figure 8 is a flow chart of the procedure beginning with a user defined 
atmosphere of known quantities of aerosol optical depth (6), water vapor 
densities (Pw), and surface ‘emperatures (Tsurf). Three programs are run for 
each data set. Both LOWTRAN7 and it's filter function program, LOWFIL7, 
make use of formatted files (TAPEX.DAT) for input, output, and to pass coded 
information between the programs. The program RETRIEVE was specifically 
designed for this study, to compile data in a condensed format. 

LOWTRAN7 is supplied a modeled atmosphere, sun/satellite/earth 


geometry, as well as a desired wavenumber («) band. Wavenumber dependent 
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TAPE6.DAT 
Broad Band Radiances 


TAPES.DAT 
NOAA-7 AVHRR 
Chan. 1,2 
Chan. 4,5 
Sensor Response Funcuons 








TAPES.DAT 


¢ 850 km satellite height 

¢ 10 km coposphere 

¢ Background scratosphenc exancaon 

e User defined aerosol scanernng phase funcnon 

¢ Varied column water vapor (x11) 

¢ Varied, user defined aerosol extinction 
coefficients (x10) 


¢ 260°K troposphenc temperature 

e Varied surface temperature (x4) 

¢ NOAA-7 AVHRR frequency windows (x2) 
Chan. 1,2 
Chan. 3.4 












TAPE7.DAT 
Broad Band Filter Input 


LOWFIL7A/B 


TAPE6.DAT 
Filtered Radiances: 
1) Total Integrated 
2)Filter Normalized 














RETRIEVE.DAT 
Chan. 1 Radiance 
Chan. 2 Radiance 


L1/L2 


Chan. 4 Brightness Temp 


Chan. 5 Brightness Temp 


T4-T5 = T4-T5 


Fig. 8 Flow chart of model run. 


radiances are calculated considering scattered, absorbed, and emitted radiation as 
it would be “seen” by a perfect sensor at the satellite location. Figure 9 shows 
the AVHRR sensor response function curves (P(k)) that LOWFIL7 combines 
with coded information derived from the previous LOWTRAN7 run to 
calculate the effective filtered radiances as a function of wavenumber or 


wavelength. The sensor weighted total integrated radiance 1s then determined. 


Ltot = [L(x)@(«)dx (3.1) 


LOWFIL7 was modified for use in this study to calculate the integral of the 
sensor's response function which is then used to calculate a normalized radiance 


(LOWFIL7A/B).2 


JLeK)@(K)dx 


(3.2) 
Jo(«)dx 


Li(chan) = 


The program RETRIEVE reads the output files from LOWFIL7A/B and 
writes only the normalized radiances for channels 1, 2, 4, and 5 (L1,L2,L4,L5). 
The ratio of channel 1 to channel 2 radiances is calculated to determine the 


particle size index parameter (S12) (Frost 1988). 


S1i2=7..=1. (3.3) 


2LOWFIL7 normalizes by , Pmax instead of Jocxdx. 


3Unlike Frost, this study does not first remove Lr, 
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Fig. 9 NOAA-7 AVHRR filter function response curves. 
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Channel 4 and channel 5 radiances are converted to brightness temperatures 


(Tp). 


C2K 
Tp(4,5) = — ee (3.4) 


Nene + L(chan 4,5) 


Where Tp4,5) 1S the brightness Temperature (K) for the radiance L(4,5) 
(w/cm2-steradian-cm-!), x is the central wavenumber of the channel filter 
(Kchan4 = 927.22 cm-! and Kchans = 840.872 cm-!), and Cy and C2 are constants 
(Cy= 1.1910659x10-12 w/cm2-steradian-cm-4 and C2= 1.438833 cm-! K) 
(NOAA Polar Orbiter Data 1986). Finally, the brightness temperature 
difference between Tpq4 and Tps is calculated (T4-TS). 


C. INPUT CHARACTERISTICS 

The modeled atmospheres were restricted to cloud-free marine environments 
where backscattered radiance is dominated by aerosol particles (Durkee 1984). 
Since the satellite senses the effect of vertically integrated parameters, the models 
used were homogeneous from the surface to the tropopause. For ease of 
integration, the model tropopause was set to 10 kilometers (the average height of 
the mid-latitude tropopause is about 11 kilometers (Stull 1988)). From the 
tropopause to the satellite height, only background extinction due to meteoric 
dust is represented from 30 to 100 kilometers. The following are the important 


user-defined parameters: 


aerosol optical depth (Oq(A))iecece eet can sss varied 
column water vapor (W)........<ceere rere ee eee ie ons eiowesvasiads varied 
tropospheric airgtemperaturE aT air) ee --seeee-e «>: <--.--osssereeees--- Bs 260 K 
surface temperature alist rf) acts cc-2 ee oo ceisoee ereeineeo aes varied 
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SmACe AlDEMON (SAMI). ..2<..:c2000.cs0ssssesscecessecnstescoseedds sasderses 0.0 (blackbody) 


SELSHINMS “eid 1a (() 9 Ei ape eine eas nents 7a ne Fen nee ee Sa 850 km 
Seo Go BL eer rere ate nna ot Gece Manan pea ac onde Mebtne masaa baa healers spring-summer 
marine phase function (P(@))..........0.0000..000.. from Shettle and Fenn (1979) 
wavenumber windows, resolution 
O11 OP ee ee ee 8400--20050 cm-!, 50 cm-! 
Cran AGS aaaetesecctennetare eee 770--1005 em-!, 5 cm-! 


Molecular constituents such as 03, CH4, N20, CO, CO2, 02, NO, SO2, NO2, 
NH3, and HNO3, used LOWTRAN7 model default profiles for midlatitude 
Summer and were kept constant for all model runs. The sun/satellite geometry 
remained the same throughout all the model runs with the satellite and sun in a 
straight line perpendicular to the earth's surface. 

1. Optical Depth Profiles 

Satellite windows in the red-visible and near-infrared sense the result of 
vertically integrated parameters. The most significant influence on retrieved 
radiances is the vertical integral of Oext or 6. To study this effect, ten profiles 
of 5 with constant Oext to 10 km were designed. 


Using the relationship 


and an Oscat of .05 km-! at A = .55 yum, a profile of Gscat as a function of A was 
derived (Oscat(A)). After integrating over the 10 km height to determine 6, 
Figure 10 shows this curve (M2) and four multiples of this profile that complete 


the first set of optical depth profiles. 
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Fig. 10 Multiple optical depth profiles (M). 
Figure 11 shows a second set of profiles created using the same relationship and 
adding equal quantities to create parallel profiles. 
2. Water Vapor 
All the water vapor was contained in the homogeneous troposphere (0 
to 10 km). The maximum value was selected from the tropical atmosphere 
model included in LOWTRAN7 (3.322 gmcm-2). User-defined profiles are 


input as water vapor mass density (pw) (gm m-3), which after integrating in the 
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Fig. 11 Parallel optical depth profiles (P). 
vertical gives column water vapor (w) of equivalent magnitudes in units of gm 


cm-2, 


H 
wW= jew dz (3.6) 


Ten intermediate quantities between 0 and 3.322 gm cm were selected with 0 


included as the eleventh. 


3. Surface Temperature 

T4-TS5 is a measure of water vapor (Jacobowitz and Gruber 1990) and 
will be compared with the user input water vapor concentrations, as an empirical 
method of estimating w. Ts5 is more sensitive to water vapor effects on 
atmospheric transmittance while T4 senses more of the near-blackbody, long- 
wave radiation from the surface. In order to demonstrate relationships, the 
model required a warmer surface temperature than the atmosphere when the 
troposphere is isothermal. With equal atmospheric temperature (Tatm) and 
surface temperature (Tsyrf), the measured brightness temperatures would be 
identical. Also, surface temperatures colder than the temperature at 10 km are 
unrealistic. Therefore, while the tropospheric temperature was held constant, 
the surface temperature was varied to achieve four representative modeled 


Tatm/T su rf relations hips: 


Tati ....cecccceeeee 260°K 
Tourf(W).....0000 270°K 
Tsurf(X)........... 280°K 
Tsurf(y).....c000e 290°K 
Tsurf(Z)........0. 300°K 


D. PROCEDURAL SUMMARY 

The model was run at the Naval Postgraduate School's (NPS) Interactive 
Environmental Digital Analysis Laboratory (IDEA LAB). A complete set of 
input files consisted of two sets of five 6 profiles, with eleven w variations, four 
surface temperature conditions, and two wavenumber windows. These 880 files 
(2 X 5 X 11 X 4 X 2) generated 1760 intermediate and final output files, with 


the significant data ‘retrieved’ and collated in one resultant array. All the 
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calculations were done as a large batch job requiring about 4.5 hours to run. 
The job was completed in both the single and multiple aerosol scattering modes 


of LOWTRAN7. 


pL, 


IV. RESULTS 


A. THE EFFECT OF AEROSOL OPTICAL DEPTH VARIATIONS 
ON S12 


In general, an increase in wavelength dependant aerosol optical depth 
(da(A)) resulted in an increase in Ly and L2 while S12 values decreased. Since 
S12 1s generally greater than one (Ly>L2). Increasing the radiance in both 
sensor windows by similar magnitudes has a greater influence on the smaller L, 
driving S12 toward one. This relationship will remain true because of the 
Rayleigh contribution, which is greater for L1, 

1. Multiple 5a(A) Scheme 

Increasing 5g profiles result in decreasing S12 ratios, with the rate of 
decrease diminishing due to the effect of increasing dq(A) slopes. Figure 12 
shows the S12 variations resulting from the multiple dq profiles of Fig. 10. 

2. Parallel 5a(A) Scheme 

By constraining the curves to be parallel, an increase in dq is an equal 
increase at all wavelengths. The S12 ratios decrease at a more rapid rate than in 
the previous scheme since a change in slope is no longer limiting the decrease in 
S12. Figure 13 shows the S12 variations resulting from the parallel 5g profiles 


of Fig. 11. 


B. THE EFFECT OF SURFACE TEMPERATURE VARIATIONS ON 
THE T4-T5 REPRESENTATION OF W 


Although the modeling accomplished in this study was not designed to prove 
Dalu's technique for determining vertically integrated water vapor density (w) 


from T4-TS, some significant aspects of his theory are supported by the results. 
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Fig. 12 S ,2 versus column water vapor (w) for the multiple optical 
depth scheme (M). 


Figure 14 shows the T4-Ts5 values derived from modeled atmospheres of known 
w given a 260K isothermal air temperature (Tair) and various Tsurf conditions. 
The slopes of the T4-Ts varied from Dalu's estimate, yet the linear relationship 
of Eqn. 2.13 seemed evident. 

The variation in the slopes is a result of the assumptions made in the model 
atmospheres. The use of an isothermal 10 km atmosphere is reasonable for 


NOAA-7 AVHRR channel | and 2 retrievals since these wavelengths 
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Fig. 13 S12 versus column water vapor (w) for the parallel optical 
depth scheme (P). 


(red-vis, nir) sense the vertically integrated aerosol backscattering with little 
height dependance. Since absorption plays a big role in channel 4 and 5 
brightness temperature retrievals in the actual atmosphere, Tp4 and Tp§ results 
are significantly influenced by vertical temperature variations. The isothermal 
assumption used in this model study, therefore, causes quantitative T4-TS5 results 
to be unrepresentative of the real atmosphere. However, qualitative results such 
as the linear relationship are retained. Table 2 compares the slopes derived from 


applying Dalu's theory to each Tair/Tsurf variation. 
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Fig. 14 T4-Ts5 versus column water vapor (w) for a variety of 
surface temperatures. 
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TABLE 2 


MODEL OUTPUT COEFFICIENTS APPLIED TO DALU'S 
THEORY (w = A(Ta-T5)) 


Tair/Tsurf (K/K) A (g K-! m'}) 
260/300 7615 
260/290 9944 
260/280 1.4995 
260/270 3.6653 

Dalu 1.9600 


C. THE W CORRECTION TO S12 
In order to present the correction to satellite retrieved S12 values required 
due to the presence of water vapor, the results will be displayed first graphically, 
then in a new parameterization that can be applied to the aerosol scattering phase 
functions used in retrieval of aerosol characteristics from satellites. 
1. Graphical Application 
RETRIEVE lists the following results for each modeled atmosphere: 


n Parameters: Variables: 
Li 
Sa(A) 7 =>S12 
Tsurf . 
bi ng =>T4-T§ 


A set of output variables was calculated for each variation of input parameters. 


The combined results are displayed in Fig. 15. 
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To use this figure to determine S12 corrected for w: Begin with NOAA-7 AVHRR retmeved data (S12, 
T4-Ts). Read T4-T§ value on lower half of plot and look for the intersecuon with Dalu curve. This 
will determine column water vapor. Move up (parallel to the grid lines), to the satellite measured $12 
value. Follow the contours of the nearest curve to the Sj? axis and read S}2(dry). This will be the S}2 
with the effects of water vapor in the column removed. 
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Fig. 15 S 12 correction for column water vapor. 
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Both the multiple and parallel S12 results are displayed and seem to follow a 
similar relationship to w. Note that M2 and P1 curves are coincident. They use 
the same input da(A) values that formed the basis for each set of curves. 
2. The Reformulated Parameterization 
Looking at the S42 values plotted for varying w and da(A) curves, 

there appears to be a functional relationship. It is apparent that the water vapor 
influence of S42 is not dependent on the shape of the optical depth profile, but 
was a function of w and the magnitude of $12 when w = 0 (S42(dry)). Figure 
16 shows that when the ratio of S42:S42(dry) was plotted as a function of w, the 
shape of the curve was conserved for all model outputs (multiple and parallel). 
The slight variation at higher w is partly due to computational round off errors. 


The curve that best fits the data is: 


S12(dry) = (1 + .0332 Vw) S12 (4.1) 
with r2 = 0.995. 


From Eqn. 4.1, measured $12 and estimated w from T4-T5§ determines 


S12(dry). 


D. SINGLE VERSUS MULTIPLE SCATTERING 

All data sets were run with both single scatter and multiple scatter 
approximations. The results were calculated to the seventh decimal place and 
were precisely the same for either case. The model, therefore, supports Frost's 
use of a simplified form of the Radiative Transfer Equation (RTE) that assumes 


single scatter (Eqn. 2.1). 
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Fig. 16 S12 normalized by S12(dry). 
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V. CONCLUSIONS 


Present techniques for satellite retrieval of aerosol characteristics using the 
AVHRR carried on the NOAA series platforms require that analyzed pixels be 
free of clouds. This is reasonable since liquid water effects on radiative 
processes are quite variable and difficult to measure. Clouds at any level in a 
satellite scene can obscure the sources it was designed to sense, or saturate the 
sensor with reflected radiance. It is understood that water in a gaseous form 
would still be present, but its influence on L in the wavebands most effectively 
used to estimate aerosol parameters (Chan. 1 and Chan. 2) has been assumed to 
be insignificant in previous studies. 

Although aerosol scattering is still the dominant influence on L in these 
wavebands, this study has shown that water vapor does introduce an error, and 
that the error can be removed. This does not invalidate the present techniques, 
but provides a significant improvement that can be applied through the use of the 
variable scattering phase function introduced by Frost (1988). Since aerosol 
variations can be of global scale, even very small changes may have global 
ramifications. 

Through the use of the LOWTRAN7 atmospheric propagation model, this 
study has shown support of the following techniques and theories presently 


accepted: 


¢ The concepts of Frost's particle size index parameter (S12) using 
AVHRR channels 1 and 2. 


¢ Frost's use of a simplified radiative transfer equation based upon an 
assumption of single scatter. 


¢ Dalu's water vapor retrieval technique using AVHRR channels 4 and 5 
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By introducing varying quantities of water vapor to the modeled 


atmospheres of differing optical depth (5) and Tsurf/Tair profiles, additional 


information has been generated: 


¢ The error introduced to $12 due to water vapor is demonstrated and 
quantified. 


¢ The error has been shown to be a function of $12(dry) and column 
water vapor (W). 


The influence of optical depth variations can be isolated from the variation 
due to w through the use of the variable scattering phase function (P(©)) 
proposed by Frost, (1988). Using Eqn. 2.1, $12 with 2.0gm cm-2 of water 
vapor produces a 15% error in P(©) and therefore retrieved aerosol optical 
depth (dg). 

By using a model to determine water vapor's influence on satellite-measured 
aerosol characteristics, the parameters normally estimated can be strictly 
controlled. As known quantities, their theoretical relationships to satellite- 
measured quantities can be observed. The next step is to apply the correction 
discovered here to empirical satellite data with verifying in situ data sources such 


as aircraft or ships. 
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